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We demonstrate that a model of disordered Anderson lattices can account for many non-Fermi liq- 
uid features observed in a number of Kondo alloys. Due to the exponential nature of the Kondo 
temperature scale, even moderate disorder leads to a rather broad distribution of Kondo temper- 
atures, inducing strong effective disorder seen by the conduction electrons. Spins with very small 
Kondo temperatures remain unquenched and dominate the low temperature properties. The model 
predicts logarithmic divergences in thermodynamic quantities at low temperatures. We also find 
a linear temperature dependence of the resistivity, a feature that remained a stumbling block in 
previous theoretical attempts. We argue that for realistic amounts of disorder, such marginal Fermi 
liquid behavior is a very robust feature of disordered Kondo alloys. 



PACS Nos. 71.10.Hf, 71.27. +a, 72.15.Qm 

Non-Fcrmi liquid (NFL) behavior in metals represents 
one of the key unresolved issues in condensed matter 
physics. There exists by now a large class of non- 
magnetic metallic f-electron materials which do not be- 
have as Fermi liquids at low temperatures . In some 
of them, the proximity to a T = quantum critical 
point appears to be the origin of the anomalous behav- 
ior |§|||[§- However, in several other cases, NFL behav- 
ior only occurs when the system has been sufficiently al- 
loyed so that it is not close to any phase boundary. This 
is the case of the alloys UCus-sPds p],p[, Mi_ x U x Pd3 
(M = Sc,Y) !§, Lai_a.CeBOua.aSia g, Cei_ x Th x RhSb 
H Ui_ x Th x PdaAl 3 0. In all of these systems the spe- 
cific heat varies as C(T)/T as aln(_b/T) and the resis- 
tivity is linear with a large zero temperature intercept 
p(T) w po(l — T/Tq). The magnetic susceptibility has 
been often fitted by a logarithm or a weak power law. 

Some attempts have been made to explain the anoma- 
lous low-temperature properties based on exotic one- 
impurity mechanisms, such as the quadrupolar Kondo 
model |lC(| . Inconsistencies with the predictions of the 
model for the resistivity (w \T) and in an applied mag- 
netic field in some of these systems, however, invite the 
consideration of other mechanisms for NFL behavior. 

Quite generally, the large residual resistivity of these 
systems together with their alloy nature immediately sug- 
gest that disorder could be significant. In an important 
recent study g], the strong broadening of the copper 
NMR line of UCu5_. E Pd. E (x = 1 and 1.5) has provided an 
independent indication of the essential role played by dis- 
order in at least one of these compounds. These results 
suggested the presence of strong spatial fluctuations in 
the characteristic Kondo temperature Tk of the local mo- 
ments Indeed, by using a model distribution func- 
tion P(Tk) and well-known single-impurity results, they 



were then able to quantitatively describe the low temper- 
ature thermodynamic properties (specific heat and mag- 
netic susceptibility) as well as the NMR linewidths. The 
proposed picture implicitly assumes independent local 
moments which is usually sufficient for understanding the 
thermodynamics of most heavy fermion compounds. Of 
course, in the context of transport in concentrated Kondo 
systems, such an assumption appears to be unjustified, 
since it cannot be reconciled with the well established 
coherence effects at low temperatures. 

The central question addressed in this letter is whether 
disorder effects can explain not only the thermodynam- 
ics, but also the anomalous transport in these systems. 
We will formulate a theory appropriate for concentrated 
magnetic impurities, which can describe the coherence 
effects in the clean limit. We will show that correlation 
effects will strongly enhance any extrinsic disorder, gener- 
ating an extremely broad distribution of Kondo tempera- 
tures. This leads to the destruction of coherence and, for 
sufficient disorder, to the breakdown of Fermi liquid be- 
havior. The low temperature properties can be viewed as 
resulting from a dilute gas of localized elementary excita- 
tions: those Kondo spins that remain unquenched. This 
picture of dirty Kondo lattices, similar in spirit to the 
original Landau description of simple metals, provides 
a clear theoretical underpinning for one route towards 
marginal Fermi liquid behavior. 

We start with a disordered non-dcgcncratc infinitc-U 
Anderson lattice model 
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where, Cka destroys a conduction electron with momen- 
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turn k and spin a from a broad uncorrelated band with 
dispersion e(k) and half bandwidth D and fj a destroys 
an f-electron at site j with spin a. The infinite-U con- 
straint at each f-orbital is assumed (nj < 1). The 
on-site energies E- and the hybridization matrix ele- 
ments Vj are assumed to be distributed according to 
some distribution functions P\{Ef) and Pity). In the 
Kondo limit, the local Kondo temperature is given by 
T Kj = Dexp(Et /(2p V 2 )) {p w ±) and will be corre- 
spondingly distributed. Because of the strong scattering 
off the f-sites, disorder in f-parameters is dominant and 
we will thus neglect other types of disorder in the c-band. 

To analyze the properties of our model, we focus on 
the dynamical self-consistent theory of strong correlation 
and disorder [|l2| . The problem can then be reduced to an 
ensemble of one-impurity problems in a self-consistently 
generated self-averaging bath of conduction electrons. 
The equations simplify considerably in the case of a semi- 
circular conduction density of states, where the ensemble 
of impurity problems is governed by the action 
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where the infinite-U constraint is implied and 
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uj - t 2 G c (uj) 



(2) 
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Here, t is the hopping parameter and G c (lo) is the 
disorder-averaged local conduction electron Green's func- 
tion. The latter is determined self-consistently by 
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where 



Here 
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) av denotes the average over disorder and 
£™ p (w) is the f-electron self-energy derived from the im- 
purity model of Eq. (g). In absence of disorder, these 
equations reduce to the dynamical mean-field theory of 
the Anderson lattice [H| , while for U — they are equiv- 
alent to the CPA treatment of disorder |JL4| for the con- 
duction electrons Jl5[ . In general, the theory is exact in 
the limit of large coordination. Once G c (cu) has been de- 
termined, the conduction electron self-energy S c (w) can 
be obtained from 
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where p (e) = y/l - (e/2t) 2 /nt. 

Let us analyze the qualitative behavior of <&j(u>). From 
the Fermi liquid analysis of the impurity problem, it is 



* e 0.2 




0.05 0.10 
Width of E f distribution 

FIG. 1. Scattering rate as a function of the width of the Ef 
distribution. The parameters used are shown in the figure. 
The strong correlations in the f-shell produce an enhanced 
effective disorder. 

well known that E^ p (u; = 0) is a real quantity at T = 
Jl6| . Therefore, one can write 



*i(" = 0) 



V? 



E 



f 



Re[E£"(0)] 



(7) 



^^(O) measures the scattering strength at site j at the 
Fermi level. In the clean limit, will be the same 

at every site and Eqs. (0) and (g) give S c (w) = $(ui). 
In this case, S c (o; = 0) is a real quantity, reflecting the 
coherent nature of the DC transport at zero temperature. 

By contrast, when the system is disordered, a distri- 
bution of scattering strengths $j is generated, strongly 
affecting the transport properties. By applying the large- 
N mean field theory to the impurity problems at zero 
temperature, we have solved the self-consistent problem 
defined by Eqs. (|-|). The resulting scattering rates as a 
function of the width of the Et distribution (for a fixed 
uniform value of V) are shown in Fig. [TJ. Similar results 
are obtained for a distribution of V values holding Ef 
fixed. For the residual resistivities reported for the NFL 
alloys, e. g. UCus-^Pd^, M, one can estimate Dt w 3—5. 
Due to the strong f-shell correlations, rather large scatter- 
ing rates can be generated by a small disorder strength in 
f-parameters (see Fig. [3] and Eq. ^). Comparable amounts 
of disorder, in the absence of correlations, cannot produce 
these large resistivities. 

Thus, with sufficient disorder, scattering off the f-sites 
becomes incoherent and the resistivity assumes a monot- 
ically decreasing dependence on temperature, resembling 
the single-impurity result. The actual scattering rate, 
however, requires the solution of the full set of Eqs. (S- 
§). If P(T K ) is broad enough, the low-temperature de- 
pendence can be non-trivial. To analyze that, it is useful 
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to rewrite the above equations in terms of the impurity 
T-matrix Tj mp (w) = V?Gfj(w), where G fj (w) is the f- 
Green's function computed from the action in Eq. (||). 
We find 



G e (tj) = 
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and, from Eq. 
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We now raise the temperature slightly from to T and 
denote corresponding variations by St- Then 
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where 
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Here, the temperature dependence of the self-energy is 
expressed in terms of the temperature dependence of 
the disorder-averaged T-matrix. In general, the self- 
consistency condition couples different frequencies, as 
seen in the integral equation (10b). However, the leading 
low temperature behavior is determined only by the lo = 
component of the averaged T-matrix, so in the following 
we concentrate on this object. 

Fig. shows the result of averaging the imaginary 
part the single impurity T-matrix over the distribution 
of Kondo temperatures deduced from the experiments of 
Ref. For the single impurity dependence, we used a 
simple scaling form with the correct asymptotic behavior 
at high and low temperatures. The dependence is linear 
at low temperatures. 

It is easy to understand the origin of the linear behav- 
ior and why it does not depend on the detailed shape of 
the curve. We will focus on the imaginary part of the im- 
purity T-matrix since it gives the dominant contribution. 
It has the following scaling form 
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where So is the phase shift at T 
has the following asymptotics 
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0. The function t(x) 
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FIG. 2. Temperature dependence of the imaginary part of 
the single impurity T-matrix averaged over the disorder dis- 
tribution appropriate for UCus.sPdi.s, as determined experi- 
mentally in Ref. j|. The inset shows the same quantity over 
a wider temperature range. 
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X < 1 
X > 1 



where a and (3 are universal numbers. It follows that 
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Now, for a fixed temperature T and as a function of 
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F(T/T K ) 
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F(T/T K ) is strongly peaked at T K m 0, decays as 1/T| 
and has a width of order T (see Fig. ||). For low T, it 
can be written in terms of a delta function of Tk , hence 
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where a = J dxF(l/x) 
this yields 

<5t5L ~ - 



When inserted into Eq. (10a 
iaP(0) 
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Therefore, the low temperature dependence only 
probes the P(TfA distribution at low values of Tk as 
is clear from Fig. S. In that region, P(Tk) can be taken 
to be a constant and the temperature can be scaled out of 
the average, yielding the negative linear term. As long as 
the distribution of Kondo temperatures is wide enough 
so that P(0) is appreciable, there will be a sizable lin- 
ear range. For sufficiently weak disorder, P(0) is zero or 
negligible and Fermi liquid behavior is recovered. 
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FIG. 3. Comparison of the experimentally determined dis- 
tribution of Kondo temperatures of the alloys UCus-^Pda; 
(x = 1,1.5) (from Ref. |j) with the function F(T K ,T) de- 
fined in the text. The function F(Tk,T) only probes the 
Tk = value of the distributions at low T. 



Physically, it is clear what is happening. As the tem- 
perature is raised, a few diluted spins with Tk < T are 
unquenched and cease to contribute to the resistivity. 
The linear term essentially counts the number of liber- 
ated spins. Since this number is small at low tempera- 
tures, they form a dilute system of removed scatterers, 
whose effect is additive (or rather, subtractive), rendering 
our treatment of disorder essentially exact. Thus, even 
though the zero temperature resistivity is a functional of 
the whole distribution P(Tk), the low temperature linear 
behavior is a much more robust feature which depends 
only on the low Tk tail P(0). 

Within the dynamical mean field theory, it is possible 
to show that the picture of independent f-sites of Ref. ||] 
is justified for thermodyamic quantities. An argument 
similar to the ones above then gives x{T) ~ hi(Xo/T) 
and C V {T)/T ~ ln(T /T). Again in this case, the NFL 
behavior is due to the presence of very Iow-Tr- spins. 

It is important to comment on the potential limitations 
of this approach. In particular, we note that the dynami- 
cal mean field approach cannot describe the effects of the 
RKKY interactions. One could imagine that pairs of lo- 
cal moments with very low Tk could well condense into 
RKKY singlets, affecting the temperature dependence. 
However, if P{Tk) is very broad, the fraction of Iow-Tr- 
spins is very small and they will be, in general, very far 
apart, rendering the RKKY interaction less effective. 

In summary, we have investigated the effects of dis- 
order in concentrated Kondo alloys, by solving dynam- 
ical mean field equations of disordered correlated elec- 
trons. We find that disorder effects can lead to consider- 
able modifications of low temperature properties, leading 
to the breakdown of conventional Fermi liquid behavior, 
consistent with some Kondo alloys. Essentially, disorder 



fluctuations lead to the absence of a characteristic energy 
scale associated with a conventional Fermi liquid. 
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